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Executive  Summary 


This  document  is  submitted  as  the  final  report  for  EOARD  grant  073040,  entitled  "configuration  of  wireless 
cooperative/sensor  networks". 

Our  accomplishments  during  the  course  of  this  project  are  as  follows: 

•  An  optimal  power  allocation  methodology  for  single-hop  multi-branch  networks  is  developed,  so  that  the  total 
transmission  power  for  the  source  and  all  relays  is  minimized,  subject  to  the  symbol  error  rate  at  the  destination 
remaining  below  a  prescribed  value.  Simulation  results  show  that  this  method  lowers  the  overall  transmission 
power  considerably  and  results  in  significant  extension  of  the  network  lifetime. 

•  A  network  lifetime  maximization  scheme  is  developed,  such  that,  when  the  battery  power  of  a  relay  is  depleted, 
the  node  is  removed  from  the  set  of  the  cooperative  terminals.  The  signal  is  then  routed  by  utilizing  other  nodes 
are  relays  in  a  new  cooperative  formation.  Simulation  results  show  that  applying  our  proposed  technique  results 
in  a  modest  increase  in  network  lifetime. 

•  Three  configuration  algorithms  for  single-hop  multi-branch  wireless  cooperative  /  sensor  networks  are 
developed,  which  incorporate  optimal  power  allocation,  network  lifetime  maximization,  or  both.  These 
algorithms  are  referred  to  as  OP,  EP-LM  and  OP-LM,  respectively. 

•  A  method  for  optimal  power  allocation  to  network  nodes  is  developed,  such  that  the  system  error  rate  is 
minimized  subject  to  total  transmission  power  remaining  constant.  Simulation  results  for  a  network  with  5  relay 
stations  show  a  modest  gain  (around  1  dB)  when  the  power  allocated  to  each  node  is  set  optimally. 

•  An  optimal  power  allocation  methodology  is  developed  for  multi-branch  multi-hop  networks,  such  that  the  total 
transmission  power  is  minimized,  subject  to  the  SER  remaining  below  a  given  value.  The  solution  to  the  special 
case  of  single-branch  multi-hop  networks  is  also  found. 

•  SER  minimization  problem  over  allocated  power  to  each  relay  (subject  to  total  power  remaining  fixed)  is 
extended  to  multi-hop  networks. 

•  The  problem  of  selecting  the  optimal  cooperative  route  in  multi-hop  networks  is  investigated.  Here,  we  calculate 
the  link  cost  for  three  types  of  communication  links,  namely  point-to-point,  broadcast  and  multi-point  to  point 
links.  Using  these  link  costs,  the  optimal  route,  which  has  the  lowest  cost  (or  transmission  power),  can  be  found 
using  a  dynamic  programming  approach.  Simulation  results  for  a  simple  line  network  show  that  the  cooperative 
route  selection  algorithm  provides  substantial  gains  in  terms  of  the  required  transmission  power. 

•  An  alternate  routing  strategy  is  investigated,  where  at  a  given  stage,  cooperation  is  taken  to  be  only  from  the  last 
U  stages,  rather  than  from  all  the  previous  stages  (staring  from  the  source).  This  approach  results  in  significant 
additional  savings  in  terms  of  transmission  power. 
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1.  Introduction 


Space-time  codes  over  multiple  antenna  systems  provide  diversity  and  coding  gains  over  wireless  fading  channels.  When 
employing  more  than  one  antenna  at  each  node  of  a  wireless  network  is  not  feasible,  data  transmission  from  a  source  node 
to  the  destination  can  be  carried  out  using  relay  terminals.  Such  a  cooperative  strategy  in  effect  forms  a  virtual  antenna 
array  with  potentially  substantial  diversity  gains. 

In  a  cooperative  network,  the  relay  stations  could  simply  amplify  and  re-transmit  their  received  signals.  The  transmission 
power  in  such  an  amplify-and-forward  (AF)  approach  could  be  fixed  or  dynamically  adjusted  depending  on  the  message 
route.  In  contrast,  the  relay  terminals  could  detect,  re-modulate  and  re-transmit  the  data  symbols  in  what  is  referred  to  as 
the  decode-and-forward  (DF)  strategy.  Clearly,  in  this  case,  the  processing  required  at  each  node  to  detect  data  symbols 
would  require  additional  battery  power. 

Of  special  importance  in  the  development  of  cooperative  networks  is  the  design  of  efficient  and  practical  algorithms  for 
network  configuration.  This  means  that,  for  transmitting  data  from  a  source  to  a  destination,  a  set  of  terminals  need  to  be 
selected  as  relays,  so  that  a  network  QoS  measure  such  as  symbol  error  rate  (SER),  transmission  delay  or  system  lifetime 
is  optimized.  In  resource-constrained  networks  such  as  wireless  sensor  networks  (WSN),  the  advantages  of  cooperation 
can  be  further  exploited  by  optimally  allocating  the  energy  and  bandwidth  resources  among  users  based  on  the  available 
channel  state  information  (CSI)  at  each  node.  In  this  research  work,  our  goal  is  to  investigate  various  configuration 
strategies  for  wireless  cooperative  /  sensor  networks. 

In  Section  2,  we  investigate  single-hop  multi-branch  networks,  in  which  the  data  is  transmitted  from  the  source  to  relays, 
which  in  turn  retransmit  the  information  to  the  destination.  We  will  extend  the  work  that  we  had  previously  reported  in 
[Vaz05],  [Vaz07],  in  which  the  transmission  power  of  source  and  all  relays  were  taken  to  be  equal.  A  minimum  number 
of  relay  terminals  were  selected  to  yield  a  prescribed  SER.  Here,  we  will  consider  two  strategies  to  enhance  the 
performance  of  this  algorithm.  First,  we  drive  an  expression  for  the  optimal  allocation  of  transmission  power  amongst  the 
selected  relay  nodes,  so  that  the  total  network  power  for  a  given  system  SER  is  minimized.  Next,  we  model  the  gradual 
depletion  of  the  battery  power  in  each  relay,  and  modify  the  configuration  algorithm  so  that  the  relay  terminals  whose 
remaining  power  are  below  a  certain  minimum  value  are  removed  from  the  cooperative  scheme.  This  would  extend  the 
overall  lifetime  of  the  network,  since  alternate  routes  are  then  selected  for  data  transmission.  We  will  also  investigate 
optimal  power  allocation  to  relay  stations  so  that  the  overall  system  error  rate  is  minimized,  subject  to  keeping  the  total 
transmit  power  constant. 

In  Section  3,  we  investigate  multi-hop  topologies  in  cooperative/sensor  networks.  We  will  extend  the  optimal  power 
allocation  technique  to  the  case  of  multi-hop  multi -branch  structures,  where  the  total  transmit  power  is  minimized  subject 
to  a  required  error  rate  at  the  destination.  We  will  also  study  the  problem  of  minimizing  the  system  SER  in  multi-hop 
networks  over  allocated  power  to  relay  nodes,  subject  to  total  power  remaining  constant. 

In  Section  4,  we  study  optimal  routing  strategies  in  multihop  cooperative  networks.  In  devising  suitable  routing 
techniques,  we  need  to  take  various  transmission  mechnaisms  in  such  a  setting  into  account.  These  include  point-to-point 
links,  broadcast  (point-to-multipoint)  links  and  multipoint-to-point  cooperative  links.  As  in  [Kha07],  we  will  calculate  the 
required  power  (or  cost)  of  each  of  these  communication  links,  but  unlike  [Kha07]  which  uses  signal-to-noise  ratio  as  the 
QoS  measure,  we  will  compute  the  link  costs  so  that  a  prescribed  error  rate  is  maintained.  With  individual  link  costs 
determined,  the  route  which  minimizes  the  total  communication  cost  is  selected  as  the  optimal  route  using  a  dynamic 
programming  approach.  We  will  present  simulation  results  to  demonstrate  the  energy  saving  offered  by  the  proposed 
approach  over  a  non-cooperative  scheme  in  which  data  is  simply  passed  from  one  to  node  to  the  next. 


2.  Single-Hop  Wireless  Cooperative  Networks 


2.1  Equal  Power  Allocation  (EP  Algorithm) 

In  this  section,  we  briefly  review  the  optimal  configuration  algorithm  for  a  multi-branch  single-hop  network  in  which  the 
transmission  power  of  the  source  and  all  the  relay  stations  are  taken  to  be  fixed  and  equal  [Vaz05],  [Vaz07].  Figure  1 
shows  a  generic  model  of  a  multi-branch  cooperative  network  in  which  M  relay  terminals  are  to  be  selected  out  of  N 
available  sensors.  The  distance  between  every  two  nodes  in  the  network  is  specified  as  dpq ,  where  p  or  q  could  represent 

the  source,  s,  destination,  d,  or  relay  terminals.  We  also  denote  the  channel  fading  coefficient  for  each  link  as  h  ■  which 

are  considered  to  be  independent  complex  Gaussian  random  variables,  each  with  zero  mean  and  variance  Q  ,  and  remain 

unchanged  during  a  symbol  period.  In  addition,  all  noise  signals  are  assumed  to  be  additive  white  Gaussian  with  zero 
mean  and  variance  N  ■ 


Figure  1:  System  model 


To  obtain  an  error  rate  formula,  we  consider  an  amplify-and-forward  model  for  the  relay  nodes,  and  assume  that  the 
source  and  relay  terminals  transmit  their  signals  through  orthogonal  channels.  The  receiver  then  applies  maximal  ratio 
combining  (MRC)  on  the  M+ 1  independent  copies  of  transmitted  signal  to  detect  the  transmitted  symbols.  With  M  relay 
terminals,  the  average  symbol  error  rate  (SER)  at  high  signal  to  noise  ratios  (SNR)  is  given  by  [Rib04] : 

*  Tsd  «=  1  Ysr,  7r,d 


where 


AA(2(M  +  1)!) 


and  k  is  a  constant  which  depends  on  the  modulation  type.  In  addition,  y  is  defined  as  y  =  O  ss  /  N0  where  s5 

represents  the  average  transmit  power.  Using  the  path  loss  model  for  wireless  channels,  for  any  two  nodes,  p  and  q,  we 
can  write  Q.pq=Aldvpq,  where  A  is  a  constant  and  v  is  the  path  loss  exponent.  Thus  y  can  be  written  as 

yn  =  Ass  /(dvpqN0)  and  hence, 

„  _  C(M)  1  77  n  ^1 


pe(M) 


,M+ 1  -M+ 1 
k  Tsd 


where  F(M )  =  mr  and  m.  is  a  metric  which  is  related  to  the  location  of  the  relay  terminal.  It  is  defined  as: 

1=1 
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where  pt  =  dsr.  / dsd  ,  i) h  =  dr.d  /  d sd  ,  xn.  =  x,-  /dsi  and  yn.  =  yt  /  dsd  . 


(2.4) 


From  Equation  (2.3),  for  a  desired  system  SER  ,  pe  ,  M  relays  should  be  selected  so  that  the  following  condition  on  the 


product  of  their  metrics  is  satisfied: 
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The  above  equation  is  indeed  the  basis  for  the  configuration  of  the  network.  Let  us  suppose  that  the  desired  system  error 
rate  is  denoted  by  pej  •  Initially,  no  relay  station  is  considered.  The  resulting  SER  is  calculated  using  Equation  (2.3)  and 


if  the  required  error  rate  performance  is  not  achieved,  one  relay  terminal  will  be  utilized.  To  select  the  most  suitable  relay 
terminal,  the  metric  of  all  terminals  are  calculated  and  sorted  in  an  ascending  order.  The  terminal  with  the  lowest  metric 
from  the  available  set  is  then  selected.  The  system  SER  is  checked  again,  and,  if  needed,  additional  stations  will  be  added 
one  at  a  time  until  the  prescribed  SER  is  achieved.  The  flowchart  of  the  EP  algorithm  is  shown  in  Figure  2.  The  reader  is 
referred  to  [Vaz05]  and  [Vaz07]  for  a  through  description  and  analysis  of  this  algorithm. 


Figure  2:  Flow  chart  of  the  EP  algorithm 

2.2  Optimal  Power  Allocation  for  SER-Constrained  Cooperative 
Networks 

In  this  section,  we  derive  an  expression  for  the  transmission  power  of  each  relay  in  a  cooperative  formation,  so  that  the 
total  power  is  minimized,  subject  to  the  system  error  rate  requirement  at  the  destination  being  met.  In  [Che05]  and 
[Hua07],  the  received  instantaneous  SNR  at  the  destination  is  considered  as  the  primary  QoS.  However,  symbol  error  rate 
is  a  more  meaningful  metric  to  be  used  as  the  QoS  since  it  directly  measures  the  system  performance. 


We  again  consider  a  wireless  cooperative  /  sensor  network  with  one  source  node  s,  one  destination  node  d,  and  N  passive 
nodes  that  have  a  capability  of  serving  as  relays.  Here,  the  term  passive  is  used  to  imply  that  such  nodes  do  not  have  their 
own  information  to  transmit  and  they  are  only  used  as  relays  to  retransmit  the  source  node  messages.  Each  passive  node  is 
powered  by  a  battery  with  initial  energy  Ein.  We  consider  an  amplify-and-forward  model  for  the  cooperative  network, 
where  each  relay  station  re-transmits  its  received  signal  after  amplification.  Here,  without  loss  of  generality,  we  assume 
unit-variance  noise  at  the  destination  and  relays. 


Below  we  derive  an  expression  for  the  optimal  transmission  power  of  each  relay.  We  start  out  by  writing  the  system  error 
probability  (SEP)  at  high  signal  to  noise  ratios  when  M  relay  terminal  are  utilized.  This  expression  is  given  by:  [Rib04] 
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where  Ps  and  Pr  are  the  transmission  powers  at  the  source  node  and  the  rth  relay,  respectively,  and  C(M)  is  again 
defined  as 

n>-D 
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Given  the  availability  of  the  channel  state  information  (CSI),  the  optimal  power  allocation  problem  can  be  formulated  as: 
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Theorem  2.1:  The  solution  to  the  constrained  optimization  problem  given  in  Equations  (2.8)  and  (2.9),  (i.e., 
P\  is  unique. 

Proof'.  The  objective  function  in  (2.9)  is  a  linear  function  of  power  allocation  parameters  and  thus  is  a  convex  function. 
Hence,  it  is  enough  to  prove  that  with 

Df  :=  {Pr  <=(0,cc),r  =  \  f(Pv---,PM)  <0}  ,  f-Df  -»  ft ,  f  (P^,- •  •  ,PM)  is  a  convex  function. 

Similar  to  [Boy04],  it  can  be  readily  verified  that  /  (Pl ,  •  •  • ,  PM )  is  posynomial  function.  By  verifying  that  its  Hessian  is 
positive  semi-definite,  it  can  be  shown  that  the  function  is  convex  on  the  non-negative  orthant. 


The  optimal  power  allocation  strategy  is  given  in  the  proposition  below: 


Proposition  2.1 :  For  the  set  of  selected  relays  in  the  network,  the  transmission  power  for  the  /'th  relay  is  given  by: 
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(2.10) 


Proof.  The  Lagrangian  for  the  constrained  optimization  problem  in  Equations  (2.8)-(2.9)  is  given  by 
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For  nodes  r  =  1,...,  M  with  non-zero  transmission  powers,  the  Kuhn-Tucker  conditions  are  given  as  follows: 
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Using  Equations  (2.13)  and  (2.14),  we  have 
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for  r  =  1, . . . ,  M  . 
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From  Equations  (2.15)  and  (2.16),  we  can  obtain  the  Lagrange  multiplier  as 
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Since  strong  duality  condition  holds  for  convex  optimization  problems  [Boy04],  we  have 

Xf(p;,-,p'M)= o. 

If  we  assume  Lagrange  multiplier  has  a  positive  value,  then  f(Px  ,  •  •  • ,  PM  )  =  0 ,  which  yields: 
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Substituting  A  from  (2.17)  to  (2.15)  results  in  Equation  (2.10).  Note  that  using  Equation  (2.16)  in  (2.10),  we  can  readily 
verify  that  P  values  in  Equation  (2.10)  are  always  positive. 


2.2.1  Optimum  Power  Allocation  (OP  Algorithm) 

in  this  approach,  we  use  the  same  general  methodology  of  the  EP  algorithm,  except  that  the  metrics  are  now  modified  to 
take  the  difference  in  the  allocated  power  to  each  relay  into  account.  Following  an  argument  similar  to  that  given  for  EP 
algorithm  in  Section  2.1,  from  Equation  (2.6),  for  a  given  SER,  we  need  to  select  relay  stations  which  have  the  lowest 
metrics,  where  the  metric  for  the  z4h  node  in  the  network  is  now  defined  as 
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(2.18) 


In  the  OP  Algorithm,  initially,  the  transmit  power  of  all  nodes  are  set  to  the  same  value  as  that  of  the  source,  i.e.,  Ps.  We 
calculate  the  metrics  of  all  nodes  in  the  network  and  sort  them  in  an  ascending  order.  If  M  relays  are  to  be  selected,  we 
choose  M  terminals  with  the  lowest  metrics.  The  transmission  power  is  then  set  using  Equation  (2.10). 


2.3  Network  Lifetime  Maximization 

In  wireless  sensor  networks,  maximization  of  the  network  lifetime  is  paramount.  This  requires  that  the  battery  power  of 
all  relay  nodes  be  continuously  monitored.  Otherwise,  stations  which  are  suitably  located  within  the  network  are  over¬ 
utilized  and  are  quickly  depleted  of  their  battery  power.  This  in  turn,  terminates  the  effective  lifetime  of  the  network,  if 
appropriate  counter-measures  for  alternate  routing  of  the  signals  are  not  foreseen  in  the  configuration  algorithm. 


Many  of  the  existing  work  on  power  allocation  in  cooperative  and  sensor  networks  focus  on  minimizing  the  transmission 
power  to  meet  the  QoS  constraint  at  the  destination  [Hon07].  While  optimal  allocation  of  power  to  each  relay  does  indeed 
extend  the  network  lifetime,  we  expect  that  incorporating  additional  measures  for  monitoring  the  remaining  energy  at  all 
relays  could  provide  further  benefit.  In  this  section,  we  first  present  a  brief  overview  of  approaches  proposed  in  the 
literature  for  modeling  and  analysis  of  the  network  lifetime.  We  then  describe  our  approach  for  extending  the  lifetime  of 
the  network. 

2.3.1  Overview  of  Modeling  Lifetime  Constraints  in  WSN’s 

Various  techniques  proposed  in  the  literature  for  conserving  energy  in  ad  hoc  and  sensor  networks  can  be  generally 
divided  into  three  categories,  namely  topology  control,  power-aware  routing  and  sleep  management,  as  described  below: 

Topology  control:  Topology  control  methods  attempt  to  preserve  desirable  properties  of  a  wireless  network 
through  reduced  transmission  powers.  A  comprehensive  survey  on  existing  topology  control  schemes  can  be  found  in 
[Sta03].  In  the  scheme  proposed  in  [Vol99],  a  node  chooses  to  relay  through  other  nodes  only  when  less  power  is  used. 
Ramanathan  proposes  two  centralized  algorithms  to  minimize  the  maximal  power  used  per  node  while  maintaining  the 
(bi)connectivity  of  the  network  [RamOO].  Two  distributed  heuristic  methods  are  also  proposed  for  mobile  networks  in 
[RamOO],  although  they  may  not  necessarily  preserve  the  network  connectivity.  Two  algorithms  are  proposed  in  [Swe02, 
Vik03]  to  maintain  network  connectivity  while  keeping  the  transmission  power  to  a  minimum.  A  topology  called 
Localized  Delaunay  Triangulation  is  shown  to  have  a  constant  stretch  factor  with  respect  to  the  original  network  [Alz03]. 
Li  et  al.  propose  a  topology  control  scheme  which  preserves  the  network  connectivity  and  has  bounded  node  degrees 
[Nin03].  The  problem  of  maximizing  network  lifetime  under  topology  control  is  studied  in  [Cal03]. 

Power-aware  routing:  Singh  et  al.  propose  five  power-aware  routing  metrics  to  reduce  energy  consumption 
and  extend  system  lifetime  [Sin98].  The  implementation  of  a  minimum  energy  routing  protocol  is  discussed  in  [Dos02a, 
Dos02b].  An  online  power-aware  routing  scheme  is  proposed  in  [LiOl]  to  optimize  system  lifetime.  Chang  and  Tassiulas 
study  the  problem  of  maximizing  the  lifetime  of  a  network  with  known  data  rates  [ChaOO].  Chang  et  al.  formulate  the 
problem  of  choosing  routes  and  transmission  power  of  each  node  to  maximize  the  system  lifetime  as  a  linear 
programming  problem  and  discuss  two  centralized  routing  algorithms  [ChaOO].  In  [San04],  Sankar  et  al.  formulate 
maximum  lifetime  routing  as  a  maximum  concurrent  flow  problem  and  propose  a  distributed  algorithm. 

Sleep  management:  Recent  studies  show  that  significant  energy  savings  can  be  achieved  by  turning  wireless 
radios  off  when  not  in  use.  In  this  approach,  only  a  small  number  of  nodes  remain  active  to  maintain  continuous  service  of 
a  network  and  all  other  nodes  are  scheduled  to  sleep.  ASCENT  [Cer02],  SPAN  [CheOl],  AFECA  [XuOO]  and  GAF 
[XuOl]  maintain  network  connectivity  while  CCP  [Xia03]  maintains  both  network  connectivity  and  sensing  coverage. 
More  recently,  a  sleep  schedule  algorithm  is  proposed  in  [Mos05]  to  maximize  the  lifetime  of  network  clustering. 


2.3.2.  Network  Power  Monitoring 

In  this  section,  we  describe  our  proposal  for  extending  the  lifetime  of  a  cooperative  network.  In  most  of  the  previous  work 
on  this  subject,  network  lifetime  is  defined  as  the  time  when  one  or  several  users  are  depleted  of  energy  [Che05]. 
However,  this  definition  does  not  accurately  characterize  the  duration  in  which  the  network  operates  properly  in  a 
cooperative  system.  A  preferred  way  of  defining  the  lifetime  of  the  network  is  by  specifying  the  time  when  the  target  SER 
at  the  destination  cannot  be  achieved  with  a  certain  probability,  and  this  is  the  definition  that  we  will  consider  in  our 
proposed  technique. 

To  incorporate  the  residual  battery  life  of  each  relay  station  into  the  configuration  algorithm,  we  monitor  the  remaining 
energy  at  each  node  after  each  transmission,  and  if  this  value  falls  below  a  certain  minimum,  we  remove  the  node  from 


the  cooperative  scheme.  Specifically,  the  remaining  energy  at  a  given  relay  node,  r„  is  denoted  by  er  (n) .  After  the 
transmission  of  /7th  message,  the  remaining  energy  at  the  selected  relays  is  written  as 

e,.  (n)  =  e,.  (n -\)-Ps,i  =  1  .  (2.19) 

We  require  that  each  relay  station  in  a  given  configuration  satisfy  the  following  condition: 

er  (n)  >Ps,i  =  \,...,M  (2.20) 

If  this  condition  is  met,  the  next  data  message  will  be  transmitted  by  the  same  set  of  relay  nodes,  otherwise  the  depleted 
nodes  are  removed  and  a  new  set  of  relays  are  selected. 


2.3.3  Equal  Power  Allocation  with  Network  Lifetime  Maximization  (EP-LM 
Algorithm) 

In  this  algorithm,  we  incorporate  the  network  lifetime  maximization  strategy  into  the  EP  algorithm.  Thus,  the  transmit 
power  of  the  source  and  all  relays  are  taken  to  be  equal.  This  means  that  the  metric  of  each  node  is  defined  in  the  same 
way  as  the  EP  algorithm.  After  each  transmission,  the  residual  energy  at  each  node  in  the  cooperative  set  is  calculated  and 
depleted  nodes  (if  any)  are  removed.  This  would  then  require  a  new  set  of  relay  nodes  to  be  selected  from  the  remaining 
terminals  in  the  network.  The  operation  of  the  EP-LM  algorithm  is  described  below  and  depicted  in  the  flowchart  of 
Figure  3. 


EP-LM  Algorithm 

Initialization: 

R  =  1,  n  =  1, 

Pi  =  Ps,  for  i  =  1, .  . .  ,N 
Recursion  1: 

Calculate  hi,  i  =  1 . N 

Select  R  terminals  that  have  lowest  value  of  hi 
Sort  all  hr  such  that  hi  <  h2  <  .  .  .  <  hR 

if  the  calculated  error  from  (3.14)  is  less  or  equal  to  the  required  SER 

if  Ps<  Er(n)  for  all  r=  1 . R 

Recursion  2 

n  =  n  +  1 

Er(n)  =  Er(n  -  1)  -  Ps,  for  r  =  1,  .  .  .  ,R 
if  Ps  >  Er(n)  for  some  r  =  1 ,  .  .  .  ,R 
Remove  depleted  nodes 
N  =  N  -  number  of  depleted  nodes 
break 

R  =  R  +  1 
if  R>  N 

break 


2.3.4  Optimal  Power  Allocation  with  Network  Lifetime  Maximization  (OP-LM 
Algorithm) 

We  now  incorporate  both  Optimal  Power  allocation  to  the  relays  as  well  as  the  network  Lifetime  Maximization  in  the 
configuration  method,  which  we  refer  to  as  the  OP-LM  algorithm.  The  transmission  power  for  all  the  nodes  is  initialized 
to  that  of  the  source,  Ps.  The  number  of  relays  is  also  set  to  M.  The  metric  for  all  the  nodes  are  then  calculated  using  the 
modified  definition  of  Equation  (2.18),  and  M  terminals  with  the  lowest  metrics  are  selected.  The  optimal  value  of  the 
transmission  power  is  then  calculated.  After  the  transmission  of  each  symbol,  the  residual  energy  of  all  relays  are 
examined,  and  the  relays  with  depleted  battery  power  are  replaced  with  terminals  with  sufficient  energy.  This  procedure 
continues  until  the  number  of  nodes  with  sufficient  energy  which  yield  the  required  SER  is  less  than  the  number  of 
remaining  nodes  in  the  network.  The  operation  of  the  OP-LM  algorithm  is  shown  below. 


Select  M  terminals  with  the  lowest 
metrics 


Figure  3:  Flowchart  of  EP-LM  Algorithm 


OP-LM  Algorithm 


Algorithm  1: 

Maximal  Residual  Energy  strategy 

Initialization: 

R  =  N  ,  n  =  1, 

Pi  =  P  s,  for  i  =  1 ,  .  .  .  ,N 
Recursion  1: 

Calculate  gir ,  r  =  1 ,  .  .  .  ,R 
Calculate  the  optimum  values  of  P  r 
if  P  r  <  E  r  (n  )  for  all  r  =  1  ,  .  .  .  ,R 
Recursion  2 

n  =  n  +  1 

Er  (n  )  =  E  r(  n  a"’  1 )  a"’  Pr ,  for  r  =  1  ,  . .  .  ,R 
if  Pr  >  Er  (  n)  for  some  r  =  1  ,  .  .  .  ,R 
Remove  depleted  nodes 
R  =  R  a"’  number  of  depleted  nodes 
break 

if  R  a%0a  0 

break 


2.4.  Simulation  Results 


In  this  section,  we  present  simulation  results  to  assess  the  performance  of  configuration  algorithms  presented  in  the 
previous  sections.  The  data  bits  are  modulated  using  BPSK.  All  communication  links  in  the  network  are  considered  to 
exhibit  Rayleigh  flat  fading.  The  source-to-destination  link  is  assumed  to  have  a  normalized  distance  of  1.  We  assume, 
(without  any  loss  of  generality)  that  the  noise  at  relays  and  the  destination  node  has  a  mean  of  zero  and  variance  of  1 .  The 
QoS  requirement  for  SER  at  the  destination  is  assumed  to  be  10  \ 

In  Figure  4,  the  average  transmission  power  from  all  the  relays  are  compared  for  the  EP  and  OP  algorithms.  For  the  OP 
algorithm,  the  number  of  the  nodes  that  can  be  selected  as  relays  is  taken  to  be  N  =  14.  In  addition,  we  assume  that  the 
coefficients  Q  fall  randomly  within  a  range  of  1  to  2  with  equal  probability.  For  the  path  loss  exponent  u  =  2,  this  is 

equivalent  to  1  <  dpq2  <  Vi.  It  can  be  seen  from  the  plot  that  the  optimum  power  allocation  scheme  significantly  preserves 
the  power  consumption  in  the  network  given  the  required  SER  at  the  destination.  For  example,  with  source  power 
Ty=20dB,  the  EP  algorithm  requires  one  relay  node  so  that  the  SER  at  the  destination  is  10-5.  This  terminal  is  selected 
optimally  and  has  the  same  transmission  power  as  the  source  (i.e.,  20  dB).  The  OP  algorithm,  on  the  other  hand,  uses  14 
relays  and  the  sets  the  power  allocated  to  each  one  optimally.  By  doing  this,  the  total  power  used  for  signal  transmission 
is  substantially  reduced  to  about  9.5  dB.  Similarly,  for  Ps=  1 4  dB,  the  EP  algorithm  results  in  the  selection  of  two  relays, 
which  together  require  a  power  of  about  17  dB.  By  contrast,  the  OP  algorithm  requires  a  total  of  10  dB  for  all  the  14 
relays  to  achieve  the  same  SER  at  the  destination.  Furthermore,  at  very  low  source  powers,  (e.g.,  Ps  =2  dB),  using  the  EP 
algorithm,  we  are  unable  to  reach  the  required  error  rate  performance  no  matter  how  many  relays  are  used.  The  OP 
algorithm  does  indeed  reach  the  required  performance  by  optimal  setting  of  the  power  for  each  node. 


Figure  4:  Performance  comparison  of  EP  and  OP  algorithms,  in  a  network  with  N  =  14  potential  relays  and  SER  of  10  5 


In  Figure  5,  we  compare  the  average  network  lifetime  for  the  OP  and  OP-LM  algorithms.  Here,  the  average  network 
lifetime  with  respect  to  the  initial  energy  at  each  node  is  depicted  versus  the  number  of  potential  relays  in  the  network. 
The  initial  battery  energy  of  relays  is  assumed  to  be  equal,  i.e.  E,( 0)  =  E0  for  all  relay  stations.  Specifically,  we  take  E0  to 
be  an  integer  multiple  of  Ps,  i.e.  E0  =  100  f\.  As  expected,  the  network  lifetime  in  both  cases  increases  as  the  number  of 
relay  nodes  or  the  value  of  Ps  increases.  However,  the  plots  show  that  the  power-aware  OP-LM  algorithm  only  provides  a 
marginal  improvement  over  the  OP  technique.  The  simulation  results  reveal  that  the  optimal  power  allocation  plays  a 
significant  factor  in  reducing  the  total  transmission  power  in  the  network  and  thus)  extending  the  network  lifetime.  The 
strategy  of  monitoring  the  energy  of  the  relay  stations  using  the  technique  applied  in  the  previous  section  at  the  relays 
provides  only  a  marginal  enhancement.  We  will  explore  other  lifetime  maximization  schemes  in  the  future. 

Finally,  in  Figure  6,  the  network  lifetime  of  the  EP  and  EP-LM  algorithms  are  compared  for  different  number  of  relays 
and  finite  battery  power  at  all  relay  nodes.  Notice  that  as  expected,  the  network  lifetime  for  the  EP  algorithm  does  not 
change  as  the  number  relays  increases  since  it  is,  in  principle,  a  function  of  the  transmission  power  of  the  relays.  For  the 
EP-LM  algorithm,  on  the  other  hand,  network  lifetime  does  increase  as  the  number  of  relays  increases.  However,  due  to 
the  non-optimal  selection  of  transmission  power  at  each  relay,  the  rate  of  increase  in  the  network  lifetime  is  modest 
compared  to  the  OP-LM  case.  Generally,  the  plot  shows  that  network  lifetime  is  extended  by  applying  the  technique  of 
the  previous  section,  especially  as  the  number  of  relays  increases. 


Figure  5.  The  lifetime  performance  of  OP  (dashed  red  lines)  and  OP-LM  (solid  blue  lines)  algorithms. 


Figure  6.  The  lifetime  performance  of  EP  (dashed  red  lines)  and  EP-LM  (solid  blue  lines)  algorithms. 


2.5.  SER  Minimization  in  Power-Constrained  Networks 


In  this  section,  we  investigate  an  alternate  problem  formulation  in  which  the  system  symbol  error  rate  for  a  multi-branch 
wireless  sensor  network  is  minimized  over  the  power  allocated  to  all  nodes,  subject  to  a  constraint  on  the  total 
transmission  power. 


Let  us  assume  again,  without  any  loss  of  generality,  that  the  variance  of  noise  at  the  destination  and  relays  is  one.  With  M 
relay  terminals  operating  under  the  amplify-and-forward  protocol,  the  optimization  problem  described  above  can  be  stated 
as  follows: 
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r= 1 


(2.22) 


where  P  ,  r  =  I _ ,  M  correspond  to  the  transmission  power  at  the  relay  nodes  and  P  corresponds  to  the  total  power.  The 

error  rate  formula  given  in  Equation  (2.21)  is  valid  at  high  signal  to  noise  ratios  [Rib04]. 


The  optimal  power  allocation  strategy  is  given  in  the  proposition  below: 

Theorem  2.2:  For  the  set  of  selected  relays  in  the  network,  the  transmission  power  for  each  node  is  given  by 
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Proof.  The  Lagrangian  in  Equations  (2.21)  and  (2.22)  is 
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For  nodes  r  =  1  with  nonzero  transmit  powers,  the  Kuhn-Tucker  conditions  are 

ffL(P„-,P„)  =  ^-h(P„-,P„)  +  A=  0, 
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Using  Equations  (2.26)  and  (2.27),  we  have 

M  f 


2  1  C(M)  1  1 


p  = 


A  k 


M+ 1 


n 


p  Q  ,  Q  . 

i  sd  rd  ;=1 


1  1 

-  +  - 


V  Ppsi  Ppid  J 


,  for  r  =  1, . . . ,  M  . 


Since  strong  duality  condition  holds  for  convex  optimization  problems  [Boy04],  we  can  write 

Xf{p;,-,p*M)= o, 
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where  f  {P*  ,■  ■  ■ ,  P*M)  =  £  Pr  —  P  .  Assuming  a  positive  value  for  the  Lagrange  multiplier,  we  get 

r=\ 

M 

f(Px  ,  *  *  * ,  )  =  0  ,  which  is  equivalent  to  ^  Pr  =  P .  Thus,  using  Equation  (2.28),  we  obtain  the  following  relation: 
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from  which  can  be  calculated  as  follows: 
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Finally,  Using  A  from  Equation  (2.3 1)  in  Equation(2.28),  we  obtain  the  allocated  power  to  each  relay  as  given  below: 
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Simplification  of  Equation  (2.32)  leads  to  Equation  (2.23). 
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It  is  clear  that  values  in  Equation  (2.23)  are  always  positive.  In  addition,  note  that  when  Qs(  =  Clrd  ,  the  allocated 
power  to  all  relay  nodes  will  be  the  same,  i.e., 


Pr=  —  ,!o,r  =  l...,M 


(2.33) 


2.5.1  Simulations  Results 

To  demonstrate  the  effectiveness  of  the  technique  presented  in  the  previous  section,  we  consider  a  network  consisting  of  5 
relays  with  the  following  node  configuration: 

*L=4’Q,d=1-  (2.34) 

r 

BPSK  is  used  for  the  modulation  of  the  data  symbols.  Notice  that  Equation  (2.23)  provides  a  set  of  relations  which  need 
to  be  satisfied  by  the  optimal  power  coefficients.  To  arrive  at  the  final  values  of  the  transmit  power  for  the  nodes,  we  use 
an  iterative  procedure.  For  the  initial  values,  we  divide  the  available  power  equally  between  the  relays. 

In  Figure  7,  the  BER  for  the  multibranch  system  is  plotted  vs.  SNR  for  the  case  in  which  the  transmit  powers  are  set 
optimally  according  to  Equation  (2.23).  The  plot  of  BER  vs.  SNR  for  the  case  of  equal-power  allocation  for  all  nodes  is 
also  shown  for  comparison.  Notice  that  at  high  SNR’s  (for  which  the  error  rate  equation  is  valid),  the  optimal  approach 
yields  a  gain  of  about  1  dB.  We  expect  similar  gains  to  be  achieved  for  other  network  configurations. 

As  mentioned  above,  the  final  values  of  the  allocated  power  to  all  nodes  are  obtained  iteratively.  Figure  8(a)  demonstrates 
the  convergence  of  power  allocated  to  the  first  relay  node.  Notice  that  after  about  8  iterations,  the  allocated  power  to  the 
nodes  becomes  stable.  Figure  8(b)  the  steady-state  value  of  the  allocated  power  to  all  nodes. 


Figure  8:  (a)  Convergence  of  the  allocated  power  to  the  first  relay;  (b)  Steady-state  value  of  the  power  allocated  to  each 

relay 


3.  Multihop  Wireless  Cooperative  Networks 


In  this  section,  we  extend  the  power  optimization  technique  derived  in  Section  2.2  for  multi-branch  single-hop  networks 
to  multi-hop  cooperative  structures.  We  will  then  obtain  the  solution  for  the  (special  case  of)  single-branch  multihop 
networks.  In  Section  3.3,  we  investigate  the  problem  of  SER  minimization  over  the  power  allocated  to  individual  nodes 
when  the  total  transmission  power  is  fixed. 


3.1.  Power  Optimization  in  Multi-Branch ,  Multihop  Relay  Networks 

In  this  section,  we  derive  the  optimal  power  allocation  to  all  nodes  in  a  multi-branch  multi-hop  network,  given  a  required 
error  rate  at  the  destination. 


Let  us  consider  that  the  network  consists  of  M  branches  from  the  source  to  the  destination,  where  each  branch  is 
composed  of  Mr  hops,  r  =  1, . . . ,  M  .  Without  any  loss  of  generality,  we  assume  unit-variance  noise  at  the  destination  and 
relays.  With  relay  terminals  operating  under  the  amplify-and-forward  protocol,  the  system  error  probability  at  high  signal 
to  noise  ratios  is  given  by  [Rib04] 
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where  P0  and  Pr  m 


are  the  transmitted  power  at  the  source  node  and  the  relay  in  the  rth  branch  and  mth  hop, 


respectively.  In  addition,  C(M)  is  defined  as 
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2(M  +  1)! 

In  addition  ,  k  is  a  constant  which  depends  on  the  type  of  modulation  used  (For  M-PSK,  k  =  2  sin2(7r/M)).  For  any  two 
nodes,  p  and  q,  Q;.  m  =  1/  d'p  q,  where  d  is  the  distance  between  nodes  p  and  q,  and  v  is  the  path-loss  exponent. 

Q0  is  the  term  corresponding  to  the  source-destination  link. 


Unlike  [Rib04]  and  [Vaz05],  in  which  uniform  power  allocation  is  considered  among  the  source  and  relays,  here  we 
minimize  the  total  transmit  power  over  the  power  allocated  to  each  node,  subject  to  the  SER  at  the  destination  being 
below  a  prescribed  value.  Given  the  knowledge  of  channel  coefficients,  the  optimal  power  allocation  problem  can  be 
formulated  as 
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Theorem  3.1 :  The  solution  to  the  constrained  optimization  problem  given  above,  i.e.,  Px  ,  •  •  • ,  PM  is  unique. 

Proof.  The  objective  function  in  Equation  (3.3)  is  a  linear  function  of  power  allocation  parameters  and  thus  is  a  convex 
function.  Hence,  it  is  enough  to  prove  that  with 
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function.  Similar  to  [Boy04],  it  can  be  readily  verified  that  f({Pr  nl)(r  m)sn)  is posynomial  function.  However,  since  its 
Hessian  can  be  shown  to  be  positive  semi-definite,  the  function  will  be  convex  in  on  the  nonnegative  orthant. 


The  optimal  power  allocation  strategy  is  shown  in  the  proposition  below: 

Proposition  3.1:  For  the  set  of  selected  relays  in  the  network,  the  optimal  transmission  power  for  the  node  in  the  r,h 
branch  and  mth  hop  is  given  by 
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Proof.  The  Lagrangian  for  the  constrained  optimization  problem  in  Equations  (3.3)  and  (3.4)  is  given  by 
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For  nodes  r  =  1, . . . ,  M  ,  III  =  1, . . . ,  M  with  nonzero  transmit  powers,  the  Kuhn-Tucker  conditions  are  given  by 
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Using  Equations  (3.8)  and  (3.9),  we  have 
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Since  strong  duality  condition  holds  for  convex  optimization  problems  [Boy04],  we  have  X  f({Pr  m)en)  =  0.  If 

we  assume  Lagrange  multiplier  has  a  positive  value,  we  obtain  the  following  relation:  f({Pr  m}(r  m)gn)  =  0 ,  which  is 
equivalent  to 
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Using  Equations  (3.10)  and  (3.11),  we  can  find  the  Lagrange  multiplier  as 
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Substituting  A  from  (3.12)  to  (3.10)  results  in  Pr  m  given  below: 
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It  can  be  shown  that  the  Pr  calculated  using  the  above  expression  is  always  positive. 


3.2.  Power  Optimization  in  SER-Constrained  Multi  hop  Links 

In  the  previous  section,  we  derived  the  optimal  power  allocation  for  the  general  case  of  a  multi-branch  multihop  network. 
In  this  section,  we  consider  the  (special)  case  of  a  single-branch  multihop  network. 

We  consider  that  the  network  consists  of  a  source  node  and  M  cooperative  nodes  that  relay  the  source  message  in  M 
successive  time  slots.  Without  any  loss  of  generality,  we  assume  unit-variance  noise  at  the  destination  and  relays.  The 
constrained  optimization  problem  can  thus  be  formulated  as  follows: 
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where  P0  and  Pm,  m  =  1 , . . . ,  M  are  the  allocated  transmission  power  to  the  source  and  the  relays,  respectively.  In 
addition,  C(l)=3/4  and  k  is  a  constant  whose  value  depends  on  the  type  of  modulation  used.  (For  BPSK 
modulation,  k  =  2  .) 


Setting  the  derivative  of  Lagrangian  in  Equations  (3.13)  and  (3.14)  to  zero,  we  get 
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Using  Equations  (3.14)  and  (3.16),  we  have 
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Using  A  from  Equation  (3.18)  in  Equation  (3.16),  we  can  find  the  allocated  power  for  the  A/h  relay  in  M- hop  wireless 
system  as  given  below: 
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3.3.  SER  Minimization  in  Power-Constrained  Multihop  Links 

In  the  previous  section,  total  transmission  power  was  minimized  over  the  allocated  power  to  individual  nodes  under  the 
constraint  that  the  system  SER  remains  below  a  given  value.  Below,  we  state  an  alternate  problem  formulation  for  a 
multihop  relay  network,  in  which  the  overall  system  error  rate  is  to  be  minimized  over  allocated  power  to  nodes,  subject 
to  the  total  transmission  power  remaining  below  a  prescribed  value.  This  constrained  optimization  problem  can  be  stated 
as  follows: 
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We  will  once  again  use  the  Lagrange  multipliers  technique  to  solve  the  above  constrained  optimization  problem.  Setting 
the  derivative  of  Lagrangian  in  Equations  (3.20)  and  (3.21)  to  zero,  we  have 
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Now,  using  Equations  (3.23)  and  (3.21),  we  can  write 
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Using  A  from  Equation  (3.25)  into  (3.23),  we  can  find  the  power  allocation  at  mlh  relay  as  given  below: 
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4.  Path  Selection  Strategies  in  Multihop  Cooperative  Networks 

4.1.  Overview 

The  amount  of  energy  required  to  establish  a  wireless  link  between  two  nodes  is  generally  proportional  to  the  distance 
between  the  nodes,  raised  to  the  path-loss  exponent.  Due  to  this  relationship,  it  is  beneficial  in  terms  of  energy  saving,  to 
relay  the  information  through  a  multihop  route.  Multihop  routing  extends  the  coverage  by  allowing  a  node  to 
communicate  with  nodes  that  would  have  otherwise  been  outside  of  its  transmission  range. 

What  makes  the  problem  of  finding  the  minimum  energy  cooperative  path  in  wireless  networks  challenging  is  the  fact 
that  an  optimal  path  could  be  a  combination  of  cooperative  transmissions,  multicast,  and  point-to-point  links  [Kha07]. 
The  problem  of  finding  the  optimal  cooperative  route  from  the  source  node  s  to  the  destination  node  d ,  formulated  in 
[Kha07],  can  be  mapped  to  a  dynamic  programming  (DP)  problem.  The  state  of  the  system  at  stage  k  is  the  set  S* ,  i.e.,  the 
set  of  nodes  that  have  reliably  received  the  information  by  the  Mi  transmission  slot.  The  initial  state  So  is  simply  {.?},  and 
the  termination  states  are  all  sets  of  nodes  that  contain  d.  The  decision  variable  at  the  M  stage,  determines  the  set  of 
nodes  that  will  be  added  to  the  reliable  set  in  the  next  transmission  slot.  The  objective  is  to  find  a  sequence  {S* }  that 
minimizes  the  total  transmition  power.  [Kha07]  refer  to  the  solution  to  this  problem  as  the  optimal  transmission  policy. 
The  optimal  transmission  policy  can  be  mapped  to  finding  the  shortest  path  in  the  state  space  of  this  dynamica  system. 
The  state  space  can  be  represented  by  a  graph  with  all  possible  states,  i.e.,  all  possible  subsets  of  nodes  in  the  network,  as 
its  nodes.  [Kha07]  refers  to  this  graph  as  the  cooperation  graph. 

The  nodes  in  the  cooperation  graph  are  connected  with  arcs  representing  the  possible  transitions  between  the  states.  As 
the  network  nodes  are  allowed  only  to  either  fully  cooperate  or  broadcast,  the  graph  has  a  special  layered  structure.  Notice 
that  the  nodes  in  the  k'b  layer  are  all  the  reliable  sets  of  size  k  +  1 .  Hence,  in  a  network  with  n  +  I  nodes,  the  cooperation 
graph  has  n+1  layers,  and  the  M  layer  has  (” )  nodes. 


This  is  illustrated  in  Figure  9  for  a  four-node  network.  The  arcs  between  the  nodes  in  adjacent  layers  correspond  to 
cooperative  links,  whereas  broadcast  links  are  shown  by  cross-layer  arcs.  The  nodes  are  marked  as  s'  and  d  for  source  and 
destination  nodes,  respectively,  and  1  and  2  for  the  two  relay  nodes.  Here,  the  link  between  the  nodes  marked  as  {s}  and 
{s,  1}  corresponds  to  a  point-to-point  link.  An  example  of  a  multicast  (or  broadcast)  link  is  the  one  between  the  nodes 
marked  {s}  and  {s,  1,  2}.  Finally,  the  link  between  {s,  1,  2}  and  {s,  1,  2,  d}  demonstrates  a  case  of  multipoint-to-point  or 
cooperative  scenario.  Some  of  the  links  are  omitted  for  clarity. 


In  [Kha07],  the  signals  arrived  at  the  destination  are  combined  and  the  resulting  SNR  is  required  to  be  greater  than  a 
minimum  value.  In  this  report,  we  follow  a  similar  strategy  except  that,  instead  of  SNR,  we  focus  on  the  system  error  rate 
at  the  destination  since  this  is  a  more  direct  measure  of  the  system  performance. 


Broadcast  Link 


Layer  0  Layer  1  Layer  2  /  Layer  3 


Terminal  Node 


Figure  9:  An  example  of  a  cooperative  graph  with  4  nodes  [Kha07] 


4.2.  Calculation  of  Links  Costs 


In  this  section,  we  compute  the  link  costs  for  successful  transmission  of  information  from  a  set  of  source  nodes 
S  =  {sl,S2,. .  .,S n}  to  a  set  of  target  nodes  T  =  {tt,t2,. .  .,tm}  .  As  in  [Kha07],  we  consider  three  different  cases, 

namely  point-to-point  links,  point-to-multipoint  broadcast  links  and  multipoint-to-point  cooperative  links.  However, 
unlike  [Kha07]  in  which  the  minimum  received  SNR  is  used  for  finding  best  cooperative  routing  path,  we  will  consider 
the  maximum  allowed  error  rate  as  the  QoS  measure.  This  metric  is  indeed  more  meaningful  in  the  design  and 
implementation  of  wireless  networks.  Below,  we  consider  each  type  of  communication  links  separately. 


4.2.1.  Point-to-Point  Links 


The  simplest  case  in  multihop  transmission  is  the  case  where  only  one  node  is  transmitting  within  a  time  slot  to  a  single 
target  node.  In  this  case,  S  =  {.Sj}  ,  and  T  =  {^} ,  i.e.,  n  =m= 1.  With  Rayleigh  fading,  M-PSK  or  M-QAM  modulations 
and  coherent  detection,  the  probability  of  error  is  given  by  [Sim05] 
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(4.1) 


In  addition,  cr  represents  the  variance  of  the  Rayleigh  fading  channel,  and  //  is  the  transmitted  signal  power  -to-noise 
ratio. 


We  assume  that  the  acceptable  bit  error  rate  for  each  link  is  less  than  BERmax  .  Using  Equation  (4.1)  and  assuming 
normalized  noise  power  at  every  receiver,  the  minimum  required  power,  and  hence,  the  point-to-point  link  cost 
ZCCsj.fj)  is  given  by 
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Notice  that  in  the  above  equation,  the  point-to-point  link  cost  is  proportional  to  l/cr“  which  corresponds  to  the 
attenuation  of  the  signal  power  in  the  wireless  channel  between  and  tx .  The  link  cost,  in  this  case,  is  therefore, 
proportional  to  the  square  of  the  distance  between  Sl  and  /,  under  our  propagation  model. 


4.2.2.  Point-to-Multipoint  Broadcast  Links 


In  this  case,  S  =  {.S',  }  ,  and  T  =  .  .,tm}  (i.e.,  n  =  1,  in  >  1  ).  Thus  here,  m  simultaneous  BER  constraints  at  the 

receivers  must  be  satisfied. 

Assuming  that  omnidirectional  antennas  are  used,  the  signal  transmitted  by  the  node  .S']  is  received  by  all  nodes  within  a 
radius  proportional  to  the  transmission  power.  Therefore,  a  broadcast  link  can  be  treated  as  a  set  of  point-to-point  links, 
and  the  cost  of  reaching  a  set  of  nodes  is  the  maximum  of  the  costs  for  reaching  all  the  nodes  in  the  target  set.  The  power 
required  for  the  broadcast  transmission,  denoted  by  LC(s1 ,  ,  is  thus  given  by 

LC  (jj ,  T)  =  ma  x{ZC(jj  ,tl),LC(sl,t1),...,LC(sx,  tm  )} ,  (4.3) 

where  LC(sl,tl)  is  given  by  Equation  (4.2). 


4.2.3.  Multipoint-to-Point  Cooperative  Links 


In  this  case,  S  =  {j15  J2,...,Jn},  and  T  =  {()},  (i.e.,  n  >  1  and  m  =  1).  This  case  corresponds  to  a  setting  in  which 

multiple  nodes  cooperate  to  transmit  the  same  information  to  a  single  receiver  node.  Assuming  coherent  detection  at  the 
destination  node,  signals  are  combined  using  maximal  ratio  combining  (MRC).  Detection  of  data  symbols  is  considered 
possible  if  the  BER  at  the  receiving  node  is  less  than  a  desired  value  denoted  as  BERmax  .  The  total  transmission  power 

*PT=±Pr 
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In  the  following,  we  calculate  the  BER  at  the  receiving  node  T  using  a  moment  generating  function  (MGF)  approach. 
With  Rayleigh  fading  at  each  link,  the  MGF  of  the  i-th  path  with  the  exponential  random  variable  y  is  given  by 


Mi(-s)  =  [  e-s,/ Pi(y)dy  = 
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(4.4) 
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Since  the  all  links  are  assumed  to  be  independent  of  each  other,  the  probability  of  error  becomes  [Sim05] 
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A.  Simplified  Case 

2  2 

Let  us  first  consider  the  case  where  all  transmitting  nodes  have  the  same  distance  from  the  destination,  i.e.,  <Jj  =  <J  , 
for  i  =  1, In  this  case,  we  can  rewrite  Equation  (4.5)  as 


p.-Li 

TT  J 


7T/2 


i  ,  8M  ,2 

I  H - j —  (7 

2  sin  cp  j 


d<p  =  -  f 

TT  J 


71 1  2 


sin2"  cp 


f  2  T 

.  2 

sin‘  (p-\ - 


dcp . 


J 
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We  will  use  the  following  identity  to  further  modify  the  above  equation:  [Gra96] 
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Using  the  above  identity  and  some  algebraic  manipulations,  the  probability  of  error  can  be  written  as:  [Sim05] 
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where 
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In  high  SNR  conditions,  in  which  gjU(7~  »  1 ,  Equation  (4.8)  can  be  approximated  as 
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Note  that  Equation  (4. 10)  also  implies  that  full  spatial  diversity  of  order  n  is  obtainable. 

B.  General  Case 

Now,  we  consider  the  general  case  of  transmitting  nodes  with  different  distances  from  the  destination.  In  this  case,  no 
closed-form  solution  of  Equation  (4.5)  exists.  In  [Rib04],  the  following  expression  for  the  BER  at  high  SNR’s  is  derived: 
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where 
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is  the  t-th  order  derivative  of  the  pdf  of  the  equivalent  channel.  The  derivatives  of  p,(/)  up  to  order 


( t  -  1)  are  zero.  Equation  (4.1 1)  is  modified  using  the  following  proposition  [Rib05]  : 

Proposition  4.1 :  Consider  a  finite  set  of  nonnegative  random  variables  \y]  ,y2,...,  Yn  }  whose  pdfs  p1,p2,...,  pn 
have  nonzero  values  at  zero,  and  denote  these  values  as  p{  (0),  p2  (0), . . . ,  pn  (0) .  If  the  random  variable  y  is  the  sum 

n 

of  the  components  of  the  set  {y[,y2,...,yn},  y  =  ^  y.  ,  then  all  the  derivatives  of  py  evaluated  at  zero  up  to  order 
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(n  —  2)  are  zero,  while  the  (n  —  1)  th  order  derivative  is  given  by 

dn^p„ 


dy 


n- 1 


(°) =n^(°) 


(4.12) 


i=\ 


The  poof  of  the  above  proposition  is  addressed  in  [Rib05]. 
Using  the  above  proposition  in  Equation  (4.1 1),  we  can  write: 
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An  exponential  distribution  with  mean  has  a  value  of  1/  cr~  at  zero.  This  will  yield  the  following  approximation: 


fl(2/-l) 


n 


i 


2  gnn\ 


(4.14) 


Note  that  if  <7i  =  <T"  ,  for  i  =  the  expression  derived  in  Equation(4. 14)  is  equivalent  to  Equation  (4.8). 


To  obtain  the  link  cost  for  the  multipoint-to-point  case,  we  need  to  consider  the  following  constrained  optimization 
problem: 
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Similar  to  the  previous  cases,  it  can  be  proved  that  the  solution  to  the  above  optimum  power  allocation  problem,  is 
unique. 

Proposition  4.2 :  For  the  set  of  n  selected  relays  in  the  network,  the  optimum  transmission  powers  satisfy  the  following 
equations: 
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Proof.  The  Lagrangian  in  Equations  (4.17)  and  (4. 18)  is  given  by 

L(Pl.-;P,)=t,P!+lf(P„-,Pl:). 
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For  nodes  /  =  1  with  nonzero  transmitter  powers,  the  Kuhn-Tucker  conditions  are 
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Using  Equations  (4.21)  and  (4.22),  we  get 

P-  =zlT(n)^TT-iT,  for  i  =  1,..., n  . 
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Since  strong  duality  condition  [Boy04]  holds  for  convex  optimization  problems,  we  can  write  A  f{P{  ,  •  •  • ,  Pn  )  =  0 .  If 
we  assume  Lagrange  multiplier  has  a  positive  value,  we  have  f  (Pl  ,  •  •  • ,  Pn  )  =  0  ,  which  is  equivalent  to 
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Using  Equations  (4.23)  and  (4.24),  we  can  find  the  Lagrange  multiplier  as 
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Substituting  A  from  Equation  (4.25)  into  (4.23)  yields  the  desired  results,  i.e., 
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Theorem  4.1 :  The  values  for  the  optimal  allocated  powers, 

Equations  (4.15)  and  (4.16)  are  identical  and  are  given  by 

^(»)fr  1 
v  BER  ix  J  ' 

Proof.  The  constrained  optimization  problem  of  Equations  (4.15)  and  (4.16)  has  a  unique  solution.  It  is  readily  verified 
that  Px  =  P1  =  ■  •  •  =  Pn  is  a  solution  for  the  set  of  the  Equations  in  (4. 1 7),  which  in  turn  leads  to  Equation  (4.27). 
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The  resulting  cooperative  link  cost  LC(S,tl),  defined  as  the  total  transmission  power  allocated  optimally,  is,  therefore, 
given  by 
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4.3.  Optimal  Cooperative  Route  Selection 

The  cost  for  different  links  calculated  in  the  previous  section  can  now  be  used  for  the  selection  of  the  cooperative  route 
from  the  source  to  the  destination  as  outlined  in  Section  4.1.  As  stated  earlier,  the  problem  of  finding  the  optimal  can  be 
mapped  to  a  dynamic  programming  (DP)  problem  [Kha07]. 

Recall  that  in  a  cooperative  graph,  the  state  of  the  system  at  stage  k  is  the  set  Sk  defined  to  be  the  set  of  nodes  that  have 
reliably  received  the  information  by  the  /cth  transmission  slot.  The  set  of  nodes  that  will  be  added  to  the  reliable  set  in  the 
next  transmission  slot,  denoted  by  U k  ,  has  the  following  relationship  with  Sk+l : 

Sk+t  =Sk{JUk,  k  =  1,2,....  (4.29) 

The  objective  is  to  find  a  sequence  {Ukj  or  {.S'/,}  that  minimizes  the  total  transmitted  power  PT  given  by 

P,  =Y,LC(.Si,U„)=Y,LC(.Sl,Sl_l-Sl).  (4.30) 
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The  result  will  be  what  is  referred  to  as  the  optimal  transmission  policy  [Kha07]. 


4.4.  Analytical  and  Numerical  Results 

We  will  use  the  following  metric  for  measuring  the  energy  saving  offered  by  the  cooperative  route  selection  technique 
outlined  in  the  previous  section  over  a  non-cooperative  scheme: 


(4.31) 
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Energy  Savings  = - , 

( noncoop ) 

where  PT  denotes  the  total  transmission  power  in  each  case. 

In  this  section,  we  present  analytical  and  numerical  results  for  achievable  energy  savings  in  a  line  network.  In  particular, 
we  consider  a  regular  line  topology  where  the  nodes  are  located  at  equal  (unity)  distances  from  each  other  on  a  straight 
line  as  shown  in  Figure  10. 


Figure  10:  A  line  network 


In  this  case,  the  optimal  non-cooperative  routing  strategy  is  to  always  send  the  information  to  the  next  nearest  node  in  the 
direction  of  the  destination.  N  relays  stations  are  denoted  by  RX,R2,...,RN .  Each  R„,  n  =  1,2,...,  /V ,  transmits  its 

detected  symbol  Xn ,  while  R0  transmits  the  user  data  symbol  X0  =  X .  The  received  symbol  at  each  node  R,„  n  = 
1,2,  ...,N  +  l  is  now  renamed  as  y„  =  hn  Xn_x  +  z,„  where  hn  denotes  the  fading  coefficient  of  the  link  between  R„  and 
R„,\.  We  assume  that  relay  stations  transmit  their  signals  over  mutually  orthogonal  channels. 

The  BER  from  R„_ j  to  R„  is  denoted  as  P2  ,  n  .  The  BER  P'n’  at  node  R„  is  affected  by  all  previous  n-1  hops  and  can  be 
iteratively  calculated  using  the  following  recursive  formula:  P%  =  (1  -  P*_  ,  )  P^_x  n  +  (1  —  P^  ,  n  )P,l’_l ,  with  P')  =  0  . 
The  end-to-end  BER  at  the  destination  can  then  be  written  as  P^+1  =  ( 1  —  P!(  )  P!)  N+l  +  (1  —  P^  N , ,  )  P!)  . 

Given  the  assumption  regarding  equally  spaced  nodes,  we  require  identical  error  rate  performance,  BERmax ,  at  all  links. 
With  BPSK  modulation  for  which  g  =  2  and  C  =  1 ,  the  required  power  at  each  link  becomes: 
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and  thus  the  total  power  for  the  non-cooperative  scheme  is  given  by 
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Let  us  now  consider  the  total  power  due  to  the  cooperative  signal  transmission.  At  each  transmission  slot,  all  nodes  that 
have  received  the  information  cooperate  to  send  the  information  to  the  next  node  along  the  same  route  as  the  non- 
cooperative  path.  After  the  mtli  transmission  slot,  the  reliable  set  is  Sm  =  {.S',1, . . . ,  in)  ,  and  these  nodes  cooperate  to 
send  the  information  to  the  next  node,  (m  +  1).  The  corresponding  link  cost  (from  Equation  (4.28))  is  given  by 
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where  BERmax  is  the  required  BER  at  each  stage  for  a  reliable  data  transmission,  and  xY(n)  was  defined  in  Equation 
(4.18). 


Using  VF( m  +  1)  from  Equation  (4.18)  in  the  above  expression,  we  can  write: 
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The  total  transmission  power  for  the  cooperative  strategy  is  then  given  by 

N- 1 

PT(coop)  =  YJLC(Sm,m  + 1)  .  (4.36) 

m= 0 

We  can  now  use  Equations  (4.33)  and  (4.35)  to  obtain  the  Energy  Saving  offered  by  the  cooperative  scheme  over  the  non- 
cooperative  approach,  as  shown  below: 
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Simplifying  Equation  (4.37),  we  get 
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4.4.1.  U-Cooperation 

In  calculating  the  Energy  Saving  in  Equation  (4.38),  we  had  assumed  that  all  previous  nodes  cooperate  in  the  (m+l)th 
step  of  transmission.  We  can  reduce  the  total  transmission  power  by  considering  the  cooperation  at  a  given  stage  to  be 
from  the  last  U  stages.  We  will  refer  to  this  method  as  the  U-cooperation  scheme.  The  Energy  Saving  for  this  modified 
cooperative  scheme  is  given  by 
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In  Figure  11,  the  Energy  Saving  of  the  cooperative  network  over  the  non-cooperative  signal  transmission  is  plotted  for  the 
number  of  nodes  ranging  from  2  to  40,  and  for  the  BER  values  of  10"3  and  10"4.  The  plot  shows  that  the  significant  saving 
in  energy  is  achieved  in  the  cooperative  multi-hop  scheme,  especially  for  the  lower  error  rate  case.  The  plots  also  reveal 
that  maximum  saving  is  obtained  when  the  number  of  nodes  is  about  4  or  5.  This  indicates  that  the  U-cooperative 
approach  with  U  in  the  range  of  4-5  could  provide  substantial  additional  saving.  In  Figure  12,  the  energy  saving  due  to  a 


Energy  Saving  Energy  Saving 


U-cooperative  scheme  is  plotted,  where  U  =  4  and  U  =  5  corresponding  toBERmax  equal  to  10'3  and  10"4,  respectively. 
Notice  that,  in  both  cases,  the  energy  saving  is  about  90%  and  approaches  unity  as  the  number  of  nodes  increases. 
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Figure  11:  Energy  saving  vs.  the  number  of  nodes 


Figure  12:  Energy  saving  vs.  the  number  of  nodes  for  the  U-cooperative  scheme 


5.  Conclusions 


In  this  report,  we  proposed  techniques  for  setting  the  optimal  transmission  power  for  relay  nodes  in  single-hop  as  well  as 
multi-hop  wireless  cooperative  /  sensor  networks.  In  each  case,  this  was  done  by  minimizing  the  total  power  subject  to 
system  SER  remaining  below  an  acceptable  value.  Alternatively,  we  explored  cases  in  which  the  SER  was  minimized 
subject  to  the  total  transmission  power  being  constant.  We  developed  configuration  algorithms  for  single-hop  multi¬ 
networks  which  incorporate  the  optimal  power  allocation  methodology. 

For  single-hop  networks,  we  studied  a  network  lifetime  maximization  scheme,  such  that,  when  the  battery  power  of  a 
relay  is  depleted,  the  node  is  removed  from  the  set  of  the  cooperative  terminals.  The  signal  is  then  routed  by  utilizing 
other  nodes  are  relays  in  a  new  cooperative  formation. 

We  also  investigated  the  problem  of  selecting  the  optimal  cooperative  route  in  multi-hop  networks.  Here,  we  calculated 
the  link  cost  for  three  types  of  communication  links,  namely  point-to-point,  broadcast  and  multi-point  to  point  links. 
Using  these  link  costs,  the  optimal  route,  which  has  the  lowest  cost  or  transmission  power,  can  be  found  using  a  dynamic 
programming  approach.  We  also  studied  an  alternate  routing  strategy  such  that,  at  a  given  stage,  cooperation  was 
considered  to  only  arise  from  the  last  U  stages.  This  so-called  U-cooperation  strategy  resulted  in  significant  increases  in 
energy  saving  in  the  cooperative  case. 
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